TMU-NT940101 
January 1994 



Meson Structure in Deep Inelastic Scattering 



Takayuki Shigetani, Katsuhiko Suzuki Q and Hiroshi Toki [| 

Department of Physics, Tokyo Metropolitan University 
Hachiohji, Tokyo 192-03, Japan 



Abstract 

We study the deep inelastic structure functions of mesons within the Nambu and Jona- 
Lasinio model. We calculate the valence quark distributions in vr, K, and p mesons at the 
low energy model scale, which are evoluted to the experimental momentum scale in terms 
of the Altarelli-Parisi equation. The resulting distribution functions show reasonable 
agreements with experiment. We also discuss the semi-inclusive lepton nucleon scattering 
process with a slow nucleon in coincidence in the final state, which reveals the off-shell 
structure of the pion. 
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1 Introduction 

The deep inelastic scattering (DIS) of hadrons is one of the most powerful tool to in- 
vestigate their internal quark structure. In the last decade, the nucleon structure function 
has been extensively studied, and provided us with detailed information of the nucleon 
structure [T[|. The observed scaling violations of the structure functions are consistent 
with the perturbative QCD predictions, though at the present QCD can not predict the 
structure function itself. 

On the other hand, it is of great interest to clarify a connection between DIS infor- 
mation and the low energy quark models0. At the experimental large momentum scale, 
the virtual photon sees the hadron as a complicated object, which consists of valence 
quarks, sea quarks, and gluons. As the momentum scale becomes smaller, sea quarks and 
gluons are absorbed into valence quarks, and their degrees of freedom are substituted by 
the 'constituent quarks', whose dynamics is subject to the low energy QCD. If we relate 
the DIS data with the quark models, we can learn how the non-perturbative aspects of 
QCD, e.g. confinement and chiral symmetry breaking, reflect the behavior of the quark 
distributions at the DIS scale. Recently, such theoretical studies for the nucleon struc- 
ture function are made in terms of several effective quark models, which are supposed to 
work at the low energy scale |P, |^. Those works are based on the assumption that the 
structure functions at the low energy model scale = are obtained by calculating 
the twist-2 matrix elements within the effective models. In the formalism of the operator 
product expansion(OPE), the n-th moment of the structure function F2 at the scale 
is expanded as, 

/ dxx^-'F2{x,Q') = Ea"(Q',Qo) < O^iQl) > (1) 

where C"(Q^,(5o) is the Wilson coefficients calculated by the perturbative QCD, and 
< (9"(Qg) > is the expectation value of the local operators which are evaluated at the 
arbitrary scale Qq. In the Bjorken limit (Q^ — > 00), only the twist-2 (r = 2) term 
survives and the higher twist terms become negligible (~ O^l/Q"^)). Hence, once we 
calculate the twist-2 operators within the quark model at the low momentum scale Qq 
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where the phenomenological model makes sense, we can get the structure function at the 
experimental scale using the QCD evolution equation. Thus, the comparison with 
experiment can be made. 

Although the available experimental data of meson structure function is much fewer 
than that of nucleon, the quark distributions in mesons are also important to study the 
quark structure of hadrons. Comparing with the nucleon case, one may extract more 
directly the information of the quark-quark interaction from the structure function, since 
we can avoid solving the complicated three body problem as in the nucleon case. In 
this paper, we shall calculate the quark distributions in mesons within the Nambu and 
Jona-Lasinio (NJL) model^]. In this model, the gluon degrees of freedom are assumed to 
be frozen into a chiral invariant effective 4-point interaction in the low energy region]^. 
The NJL model demonstrates the spontaneous breakdown of chiral symmetry and the 
emergence of the Goldstone bosons. The generalized SUiK) / NJL model reproduces the 
meson properties remarkably well, in spite of the lack of confinement |^, ^. This model 
is also applied to the chiral phase transition at finite temperature and density All 
these results indicate that the NJL model possesses the essential features of QCD. Note 
that the NJL model satisfies the Lorentz invariance and certain kinematics. This is an 
advantage of our use of the NJL model, because the quark model which does not satisfy 
the translational invariance, e.g. the MIT bag model, gives an incorrect behavior of the 
quark distributions atx = a:;B = 10, This difficulty, so called the support problem, 
may be removed by the momentum projection. However, some ambiguities exist in the 
projection methods |1T0|. 



In the recent work of present authors|ll|, the pion and kaon structure functions were 
calculated in the NJL model, and our results turned out to be consistent with experi- 
ment. In that work, however, we used an artificial momentum cutoff to reproduce the 
experimentally known behavior around x ~ 1. Here, we newly fix the momentum cutoff 
procedure to reproduce the meson properties such as masses and decay constants, and 
avoid the ambiguity. 
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We also study the semi-inclusive lepton nucleon scattering as a possible way to measure 
the quark distributions in mesons. Experimentally, the pion structure function has been 
extracted from the Drell-Yan process|[l^, |13|. On the other hand, the semi-inclusive 
processes, e.g. p + >n + i" + X, are expected to allow a new determination of the 
quark distributions in the pion|n, |T^, |T^. In such processes, the virtual photon probes 
the virtual meson clouds around the nucleon, especially pion. Hence, this experiment 
will determine the off-shell pion structure function. We can make predictions for such 
semi-inclusive scattering in terms of the calculated pion structure function. 

This paper is arranged as follows. In sec. 2, the Nambu and Jona-Lasinio model is 
introduced, and its model parameters are fixed to reproduce the meson properties. Sec. 3 
is devoted to the calculation of the hadronic tensor of mesons. The structure functions of 
the pseudoscalar channel as well as the vector channel are given in the NJL model, where 
twist-2 contributions are calculated in the Bjorken limit. These quark distributions are 
evoluted to the experimental high scale with the help of the Altarelli-Parisi equation 
in sec. 4. We compare them with data from the Drell-Yan experiment. The semi- 
inclusive processes are formulated in sec. 5. The final section is devoted to summary and 
discussions. 



2 The Nambu and Jona-Lasinio model 

The Nambu and Jona-Lasinio model0 has attracted considerable interests as a low 
energy effective theory of QCD. The NJL model incorporates the chiral symmetry, which is 
one of the most important aspects of the low energy QCD, and provides the simple picture 
of its spontaneous and explicit breakdown. Recently, this model was extensively studied 
in several subjects of hadron physics, and gave successful results^, 0. The generalized 
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SU{3) f NJL model lagrangian is given by, 



^NJL — Co + 

Co = i'{il^d^ - m)i) 

-Gvmal,i^? + {^t,l,fi^?]. (2) 

Here, if) denotes the quark field with the current mass m. ti are the S'[/(3) flavor matrices, 
normalized as tr{titj) = Sij/2, and Gs, Gy the coupling constants. 

Using (H), one gets the Dyson equation for the quark propagator in the Hartree ap- 
proximation (Fig.l (a)); 

M = m-Gs<^i^> , (3) 

where M is the dynamically generated constituent quark mass. The quark condensate 
< > is expressed as, 

<**>=-.«./p^-^. (4) 

where No is the number of colors. If the coupling constant Gs is larger than a critical 
value, the quark condensate and the constituent mass have non-zero values, which means 
spontaneous breakdown of the chiral symmetry. 

The meson properties are obtained by solving the Bethe-Salpeter (BS) equation, illus- 
trated in Fig.l (b). The inhomogeneous BS equation for q-q T- matrix is written symbol- 
ically as, 

T = /C + }CJT , (5) 

where T is the T-matrix, /C is the interaction kernel obtained from (^, and the loop 
integral of a given channel. /C is decomposed into the various products of the gamma. 
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flavor and color matrices. (See ref.p] for detail.) The loop integral J is defined by, 



where Spik-, M) is the quark propagator. 
1 



Sp{k,M) 



The quark propagator contains the constituent mass obtained by solving the gap equation 
(0). The meson mass is determined as a pole of T- matrix at = rn^^^^^. The quark- 
antiquark T-matrix is rewritten as, 

r{q) = /c + Kjr 



l-JK, 



1 

q--m-^^ 



Here, Qmqq is the meson-quark-quark coupling constant. The integrals in (|^) and (||) 
diverge due to the non-renormalizablity of the NJL model. Thus, we introduce the Fermi- 
distribution type momentum cutoff function in these integrals; 

1 



E- 



' I + exp[{kl - K^) / a] ' ' ' 

Here, k'^ is the Euclidean four momentum square, and A is identified with the typical 
scale of the chiral symmetry breaking ~ IGeV. We use a ~ O.lGeV^ to reproduce the 
meson properties. This form of the regularization function is consistent with the usual 
Euclidean sharp cutoff scheme]^]. One may understand the physical implication of the 
cutoff procedure as an approximate realization of "Asymptotic Freedom" in the NJL 
model, i.e. the interaction between two quarks with the relative momentum larger than 
A is turned off, and two particles are free in such a high momentum region. Hence, the 
structure functions obtained in the NJL model exhibit the exact Bjorken scaling in the 
deep inelastic limit |T7 . 



Solving the BS equation (|^), we can obtain the meson masses, and the decay constants. 
We take the current mass of the up and down quarks m„ = 5.5MeV as a semi-empirical 
value. The remaining parameters are fixed by the pseudoscalar meson properties. The 
parameters and the resulting physical properties are tabulated in Table 1. Using this 
model, we shall describe the deep inelastic structure of mesons. 



3 Calculation of the Hadronic Tensor 



We evaluate the structure functions of mesons in the NJL model, as done in ref. ||Tl 
The calculated distribution gives a boundary condition for the QCD perturbation at the 
model scale. The hadronic tensor W^i, is written by the structure functions Fi{x) and 
F2{x) in the scaling limit; 

+ ^iP^-^-^^^)iP^-^-^^^)^^i^) (9) 

where 

F2{x) =xJ2e^^[<l^{x) + Qii^)] , Fi{x) = ^F2{x) 

qi{x) and qi{x) are the momentum distributions of i-fiavor quark and antiquark. The 
hadronic tensor is related to the forward scattering amplitude T^,^ through the optical 
theorem [Q. 

W^, = ^ImT^, . (10) 

Thus, we calculate T^^^ in the NJL model to get the structure functions. 

We first consider the pseudoscalar meson case. We compute the forward scattering am- 
plitude in the impulse approximation, which is illustrated in Fig. 2 ("handbag diagrams"). 



Here, p is the momentum of a target meson, and q the momentum dehvered by the vir- 
tual photon. We expand the matrix elements of Fig.2 as a series of l/Q^- The leading 
0((l/(5^)°) term of the "handbag dia gram" coincides with the twist-2 contribution of the 



OPE in the Bjorken limit [0. Thus, this diagram is enough for our purpose. We note 
that, though leading terms do not depend on Q"^ explicitly, they receive the logarithmic 
QCD radiative corrections, which are incorporated by the Altarelli-Parisi equation. The 
matrix elements of Fig.2 are given by. 



f d'^k 1 

T^^ = i J (^^^^^t^/^^^^^'^^-l + *^^+ ^erm), (11) 
where 

T_ = Spik - q, Mi)gpggr+i-f5SFik - q - p, M2)gpqqTA'-fr^SF{k - q, Mi) (12) 

T_ represents the contribution with an antiquark being a spectator. T+ also has a similar 
expression, where a quark is a spectator. Here, Qpgg is the coupling constant among 
two quarks and the pseudoscalar meson obtained from (^, Q the charge operator, and 
T± = (ti ± iT2)l\pi. Ml is the constituent mass of the struck quark, and M2 the mass of 
the spectator antiquark, which are solutions of the gap equation (|^) due to the dynamical 
chiral symmetry breaking. In the case of pion, we set Mi = M2 = Mu{d)- For meson. 
Ml = Mu and M2 = Ms. 



We shall carry out the integration of (|T2|) in the Bjorken limit [p^, ITT 



-^2 



Q = —q — > 00 , u = > 00 , X = : nxed. 

nips 2mpsU 

Here, x is the so-called Bjorken x, and rrips the pseudoscalar meson mass. We introduce 
the Sudakov variables; 

kf, = zp^ + -I- , (13) 
where satisfies k ■ p = k ■ q = 0. Thus, is spacelike (k^ < 0) and is effectively 



8 



two dimensional. Calculating the traces and neglecting irrelevant terms in the Bjorken 
limitll^, we obtain, 



7-(-) 



9 (27r)' 



[Ncglqq j dzdyd^ K[ti{^^ , s) + ^t2(/i^ -5)]^ 



X 



2;z 1 , 



(14) 



where 

^2(/i^s) 

and 



1 



1 



{fi^ - + iey s - +ie 

{s + Ml - 2M1M2 - , 



{k — qY = zy + 

{k — q — pY = {z — l){y — p^) + . 



Here, we change the variable y = 2mpsv{y — 1) + zp^ [^, and yU^ and s are the invariant 
masses of the struck quark and spectator. Performing the z-integral in (0), we findO, |T3], 



8 1 c 
= 9 ^^^^^^'m j dyd^f^MfJ'^ , s) + Zt2(/i^ , s)] 



X [-g^,u + -^P^,Pu + -^—{Pfiqu + Puq^^)]■ (15) 
rripsiy ^ rupsiy 



It is easily seen from (p!5| ) that the calculated structure functions exhibit the Bjorken 
scaling[|l^]. Consider now the y-integral in the complex y-plane. The s-propagator has 
a pole at y = {k^ — M|)/(l — x) + p^ + ^^/(l ~ 2;), and the /i-propagator has a double 
pole at 1/ = (Mf — H^)/x — ze/x. For x > 1 or a; < 0, all these singularities occur on the 
same side of the real |/-axis, and ( [15|) gives a zero result |19|, |l^. However, if < a; < 1, 
the integration over y no longer vanishes. Integrating ([ISD by y and taking the imaginary 
part, we get the quark distribution by use of the optical theorem (p^O]). We change the 
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integral variable to fi^, and find the following expression for the quark distribution 
function. 



q{x) oc -g 




1 2M1M2 - {Mf + M|) + p2 



x^(/x(l - x) - xM2^ - (1 - 



(16) 



Here, 6 is the usual step function which arises from the spacelike condition of k. We 
identify (p!6|) as the valence quark distribution of the pseudoscalar meson qvai{x). In the 
case of on-shell pseudoscalar meson, we use = 77^2 ^j^g Euclidean variable 



the resulting distribution shows a correct behavior — as x — 1, since the lower 
limit of the integral fiEmini— ~/^max) — ~ ^P"^ ^00 as x ^ 1 and thus the 



the distribution function is small. This smallness is due to the spontaneous breakdown of 
the chiral symmetry. In fact, the second term disappears in the chiral limit; m„ = 771^ = 
rris = 0. This form ensures the behavior xq^aiix) oc x at small x. If the chiral symmetry 
were not spontaneously broken, the second term would be as large as the first term and 
the pionic quark distribution would behave xq{x)vai oc x^ around small x. 

We also apply the same procedure to the vector meson, though the NJL model does 
not work well in the vector channel]^. Since the vector mesons, e.g. p meson, may 
be weakly coupled quark-antiquark states, the confinement of quarks is essential for the 
description of them. Hence, it is not adequate to use the NJL model for the vector mesons 
due to the absence of confinement in this model. However, such an attempt is still worth 
mentioning for an intuitive understanding of the relation between the structure function 
and the quark interaction. 

The hadronic tensor of the spin-1 vector meson is defined as. 



2 for the integration of (|r^) with the relative momentum cutoff (|^). Note that 



integral vanishes [P^. We also note that the contribution of the second term of ( [TB| ) to 
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+ -^iP, - ^q,)iPu - ^qu)F2{x)] i-g'^^ + ^) (17) 

where a, j3 is the polarization of the vector meson, and my the vector meson mass. This 
form guarantees the current conservation, i.e. pa W^^ = PfjWl^^ = 0. Contracting with the 
meson polarization vector ea, and summing over the helicity A, we find the expression 
of the unpolarized structure function; 



W^, = Y.^c.iX,pre^iX,p)W^ 

A 



af3 



myv q^ q^ 



Corresponding to T"^ in the vector channel is given by, 



= Spik - q, M)gv,,T+rSFik - q - p, M)gv,,T_^^ S^ik - q, M) . (19) 

Here, we omit the mass difference of the struck and spectator quarks for simplicity. After 
some trivial algebras, we obtain the quark distribution in the vector meson. 

xe{p^x{l- x) - xM^ - {1- x)i/) (20) 

Here, p"^ = my. We remark that the expression for the quark distribution in the vector 
meson is quite different from that of the pseudoscalar meson. In contrast to the pion case, 
the quark distribution function of the rho meson behaves xq{x)yai x'^ at the small x, 
due to the large contribution of the third term in (PU[). 
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Unfortunately, it was shown that the mass of p meson is more than two times the 
constituent quark mass, and it should decay into two quark states, within the standard 
parameters of the NJL model[§]. Thus, we use a larger value for the constituent quark 
mass (400Me^) and Gy without changing other values to obtain a 'bound state' for p 
meson. 



4 Numerical Results 

In this section, we will show numerical results for quark distribution functions with 
the use of the parameters in Table 1. We first show in Fig. 3 the quark distribution in the 
pion (|16D at the low energy model scale. The peak of the resulting distribution appears 
at X ~ 0.6, which indicates asymmetric momentum distributions in the pion; the struck 
quark carries a larger part of the pion momentum. This is due to large binding energy of 
the valence quark in the pion. Such a behavior is also suggested by the QCD sum rule 
calculation ||2^. This result is a consequence of the highly non-perturbative structure of 



the pion 

The valence quark distribution of kaon is also interesting. We can extract the effects 
of the explicit SU (3) flavor symmetry breaking from the quark distribution of kaon. The 
u- and s-quark distributions of the kaon at the low momentum scale are presented in 
Fig. 4. The heavy strange quark carries a larger fraction of kaon momentum than the 
light up (down) quark, as expected. 

We note that this low energy scale structure function has no physical meaning at 
this scale, since "real" structure function at the low energy scale receives non-negligible 
contributions from all twist operators. The calculated results in Fig. 3 and Fig. 4 play roles 
of only the boundary condition of the structure function for the evolution. 

We take the low energy hadronic scale at Qq = (CSGeV^)^, which is used in ref. ^ 



At this scale, the running coupling constant is still small; asiQD/'K ~ 0.3. Indeed, the 
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inclusion of the second order QCD corrections gives a small change for the evolution 
from our result within 10%|pTl . 



We may understand intuitively the physical meaning of this scale Ql as compared to 
the valence quark core radius of pion {r'^)core, as noted by Brown et a/. [^. In their 
analysis of the pion electromagnetic form factor, the pion consists of the valence quark 
core part, where a valence quark and an antiquark move in the small region, and the sea 
quark cloud part which is understood as the old vector meson dominance. Their value is 
consistent to our low energy scale; 

(r')core ~ (0.35/m)2 ~ l/{0.5GeVf . 



We use the first order Altarelli-Parisi equation |^ for the evolution of valence 



distributions with Aqcd = 250MeV to compare our results with experiment. 



rf(logQ2) 2n 



Here, V{x, Q ) is the valence quark distribution, and Pqq the splitting function|23| . 

We show in Fig. 5 the pionic quark distribution at = 20GeV^ by the solid curve 
with experimental data (the dashed curve) extracted from the Drell-Yan process[^. We 
find a reasonable agreement with experiment. The second moment of the valence quark, 
which is identified with a momentum fraction carried by the valence quark, turns out to 
be, 

< xu >^= 0.22 at g2 ^ 20GeV^ , 
where < xq >= jQdxxq{x). This value is remarkably consistent with experimental 
data 0.21p3. However, the calculated distribution function is almost zero at x ~ 1, and 



different from the experimental fitpj] or the counting rule prediction This shortcoming 
comes from the cutoff procedure of the model. Around x ~ 1, the struck quark has a 
very large momentum > IGeV, and the quarks with so large momenta are excluded in 
the NJL model by the cutoff. At the moment, such a high momentum quark can not 
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exist in the hadron wave function within the low energy quark model, and we ought to 
develop a model to include the high momentum correlations consistently with the low 
energy theory. 

If we vary the low energy scale Qq within 20%, the change of the quark distribution 
is rather small (~10%). As an example, we plot in Fig. 6 the result in which we use 
Ql = 0.7GeV'^ ~ ^NjLcuf The peak position of the calculated distribution is x ~ 0.5, 
which disagrees with experiment. 

We compare the pion structure function Ft^{x) with experiment in Fig. 7. Here, we do 
not take into consideration the sea quark distributions, since we do not evaluate the sea 
quark diagrams at the model scale. It shows a good agreements for x > 0.3. In the low x 
region, the structure function is dominated by the sea quark distributions. As for the sea 
quarks, we can calculate their distribution functions at the model scale as higher order 



loop corrections in the NJL model [^. We will discuss this point later. 

The quark distributions in kaon are also evoluted to the experimental scale in Fig.8. 
The resulting distributions show a similar flavor dependence with the previous results 
obtained by the Regge theory p6|, ^ . The corresponding second moments of the valence 
quarks are given by, 

< xu >K= 0.20, < xs >K= 0.24 at ^ 20GeV^ . 

This result indicates that the heavy strange quark has a larger momentum in kaon than the 
light u(d) quark. These values are to be compared with the pion case; < xu >n= 0.22. The 
total momentum carried by the valence quarks in the kaon is < xu >k + < xs >k= 0.44, 
and is almost the same as that in the pion 2 < xu >7r= 0.43. Similar results are obtained 



by the recent QCD sum rule calculation! 28 



We also show in Fig. 9 the ratio of kaon to pion valence u-quark distributions uk/u^^ at 
= 20GeV^. The experimental values are taken from the Drell-Yan experiment ||13||, and 



the analysis of the large transverse momentum ttq production processes[^. The result is 
consistent with available experiments. Note that this ratio is sensitive to the mass differ- 
ence of the constituent quark masses in our model. In fact, if we change the constituent 
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strange quark mass, the resulting ratio becomes quite different from experiment. We plot 
two cases (M^ = M„ and Mg = 2M„) in Fig. 9 for comparison. 

Finally, we show in Fig. 10 the quark distribution in the rho meson at the model scale. 
This behavior is extremely different from the result of the pion. The resulting distribu- 
tion shows a sharp peak at x ~ 0.5 due to the weak correlation in the vector channel, 
and resembles with the result of the non-relativistic static limit q{x) ~ 6{x — 1/2). Sim- 
ilar behavior is obtained in ref. Although the NJL model in the vector channel 
has some difficult ies0, the essential feature of this distribution is expected to be model 
independent. 



5 Semi-inclusive lepton nucleon scattering 

Since a direct experiment of the lepton-pion deep inelastic scattering is beyond the 
present experimental abilities, the information of the pion structure function can be 
extracted from only the /x-on pair production process, as first suggested by Drell and 
Yan|T^. Several authors considered an alternative approach to determine the pionic 



quark distributions [0, |T^, [TBI. In the semi-inclusive lepton nucleon scattering illustrated 



in Fig. 11, virtual meson clouds mainly contribute to this process in the case of a slow 
nucleon in the final state. This contribution is expected to be dominated by the one-pion 
exchange from the study of the inclusive process p + p^n + X, in which the kinematical 
condition is similar with the semi-inclusive process mentioned above; the proton beam 
carries the incident energy ~ lOGeV and the momentum of the final state neutron is 



bOOMeV 31 



It is also important to study such virtual pion processes in view of the recent exper- 



imental observation, i.e. 'the violation of the Gottfried sum rule'p2[. The pion clouds 



make substantial contributions to the violation of the flavor symmetry in the nucleon sea 



3^ . In such a virtual pion lepton scattering, the current sees the off-shell structure 

15 



of the pion. Here, we reexamine these processes by taking into account the off mass shell 
dependence of the pion structure function, which is neglected in the previous works. 

We consider the following process; (Similar expression can be developed for the neutrino 
beam.) 

e + p-^e' + n + X. 

The calculations are performed at the target rest frame with the z-axis fixed by the 
current. Following the work of Lusignoli et ai, we define kinematical variables shown in 
Fig.lipi. 



k = £"^(1, sin-?/', 0, cos?/') for the incident electron 

k' = E'i^{l, sin^ip + 9i), 0, cos{ip + 6l)) for the outgoing electron 

Pi = M(l, 0, 0, 0) for the target nucleon 

Pf = {E,psmacos P,psmasm P,pcosa) for the outgoing nucleon (21) 

Here, M is the nucleon mass, and 6l the angle between the electron momenta, and 
tanip = E'l^smOi/ {El — E'^^ cos 6l), which is chosen to provide q = k' — k to have only 
the z-component. p is the three momentum of the final state neutron. The definition of 
the Bjorken variables are as follows; 

u = q- PJM , Q2 ^ ^ 4^^^/^ gi^2 Iq^^ ^ ^ QV2Mz/ . (22) 

For the pion, we need to define its momentum and a quark momentum fraction in the 
pion. 

P^ = P,-Pf, t = P^ = 2M{M-E), x^ = Qy2q-P^ (23) 
It is easily shown that is related with the Bjorken x as; 



X q ■ Pn , E y/u^ + p 

— = 7=r = 1 — TT H TT COS a 

Xn q-Pi M u M 



(24) 
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The requirement of the kinematics imphes the following conditions with z/ — ^ oo, 



E M 

a; < Xmax ~ 1 - < 1 , cosa> (25) 

p p 

The total (unobserved) missing mass Mx of this process is given by, 



M], = {q + P^f = Q\--l)+t. (26) 

The condition of the deep inelastic scattering also requires the sufficiently large missing 
mass. In the present study, we assume that Mx is larger than 2GeV] Mxmm = 2GeV. 
Thus, the variable is constrained by the above conditions as; 



The cross section of the electron-proton process is calculated as||15||. 



X < X/Xmax <X^< — 2 ^ 1 • (27) 

V — r + -iVlx min 



da W glxN E[ 



L 



dE'^d COS Oid^p 7rQ4 4vr ME{t-mlf 
x{2sin2l0iFf(x.,Q^^) 



+ f E^, ^'^^ 2^^] g.P. ^' ^^^^ 



where 



F-(a;^, t) = -^F2-(x^, Q\ t) (29) 

is the pion structure functions, is the pion mass, and Qt^nn = 13.5 the pion-nucleon 
coupling constant. 

For the pion exchange mechanism, we examine the elementary pion with the dipole 
ttNN form factor and the Reggeized pion. In the former case, the pion-nucleon coupling 
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becomes, 

where A = IGeV is the momentum cutoff of the nNN interaction. On the other hand, in 
the latter case, we replace the pion propagator square with the triple Regge formula |]T^]. 



1 2 /2 I +COs[7ra^(t)] b(t-m^)r V^^ " 1 " l~2a^(t) /01 N 

(t-ml)^ 2sin2[7r«,(t)] h/x~l-Myq^^ ^ ' 

Here, a.,r{t) = a'{t — m^) with a' = IGeV"'^. Exponential damping factor b is chosen 
phenomenologically to reproduce the data oip + p^n + X scattering; b = 0.56GeV~^ . 

We first show in Fig. 12 the off-shell t-dependence of the pion structure function evoluted 
to = 20GeV^, using (p^, where = t and g-nqqif)- Here, g-nqqif) is almost unchanged 
from the on-shell case. This effect is neglected in the previous works[|K|. As the pion 



momentum t increases, the peak position of the quark distribution moves toward the small 
X region, and the distribution function shows a substantial reduction for x > 0.4. Around 
X ~ 0.5, the absolute value of the distribution function at t = —0.5GeV^ is almost half 
of that for the on-shell case. This decrease is crucial to understand the following results. 

Using (|28| ) with (^) or (0), we obtain the cross section with the incident electron 
energy El = 30GeV. Here, we use the distribution functions at = 20GeV^ shown in 
Fig. 12 for simplicity, since the change of distributions due to the evolution is rather 
small (< 20%) in the momentum region above = 5GeV^, which is about the same size 
as the experimental uncertainty of the pion structure function||34]]. 

For comparison, we also present the results with other plausible forms of the pion 
structure function; the experimental fit of Sutton et and the parameterization of 

Gliick et a/.|]21]]. These are chosen to reproduce the on-shell data of the pionic distribution, 
that is, they do not depend on the pion momentum. In each figure, we show our calculation 
by the solid curve. The results with the distribution functions of Sutton et al. and Gliick 
et al. are depicted by the dashed and dotted curves, respectively. 
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Concerning the sea quark distribution, we do not calculate it at the model scale, and 
thus we can not estimate its contribution. Experimental determination also contains a 



large uncertainty ||24||. Hence, we simply use the parameterization of Gliick et a/.|^ for 
the sea quark distributions in all the cases. The sea quark contribution to the cross section 
is negligible for x > 0.3. 

We show in Fig. 13 {da/ dp) with p being the three momentum of the neutron in the final 
state. We utilize the triple Regge formula (|3TD in Fig. 13 (a), and the dipole vertex (|30|) 
in Fig. 13 (b). The results depend on the type of the ttNN vertex, especially p > 0.36*6^^, 
whereas they are almost independent of the choice of the structure function, i.e. the off- 
shell effect is not evident. Thus, we can get the most suitable form of the nNN exchange 
process from these figures by comparing with the forthcoming experimental data. 

We next discuss the x dependence of the cross section {da/dx) shown in Fig. 14 (a) and 
(b). For X > 0.3, the resulting cross section shows a clear difference between ours and 
others. The result with the off-shell structure function is much smaller than that with the 
experimental fit for x > 0.3. Note that the sea quark contribution is very small for x > 0.3. 
Hence, this difference is caused purely by the change of the valence quark distribution 
function. It is understood as a convolution effect. Roughly speaking, the cross section at 
X is expressed as an integral of the pion structure function Because of (^71) , this 

integration over Xj^ is essentially carried out from a; to 1, namely, the partons with the 
momentum fraction x^^ > x in the pion do contribute to the cross section at x. Thus, the 
cross section is more sensitive to the change of the large Xj^ behavior of the distribution 
function. We recall that the reduction of the calculated off-shell distribution function 
from the on-shell one becomes large at the large x region shown in Fig. 12. Therefore, as 
X increases, the difference of the two cross sections becomes considerably large. 

Unfortunately, all of this decrease is not originated from the off-shell effect. To see this 
point, we also present the result by the dash-dotted curve in Fig. 14, using the calculated 
pion structure function with no off-shell dependence, i.e. we neglect the t-dependence of 
the distribution function and fix t = = (140Mel^)^. It gives slightly smaller cross 
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section than those of Sutton et al. and Gliick et al. for x > 0.3. This reduction comes 
from the discrepancy between our calculation and the empirical distribution function for 
the on-shell pion, already shown in Fig.5. Because our quark distribution is different from 
the empirical one at x ~ 1, this disagreement reduces the cross section. However, our two 
results, t-dependent case and t = fixed case, also show a clear difference in Fig. 14. 
This reduction of the cross section is due to the real off-shell effect, as discussed above. 



6 Summary and Discussions 

In this paper, we have studied the deep inelastic structure function of mesons using the 
Nambu and Jona-Lasinio model as a low energy effective model of QCD. The calculations 
of the leading twist operator are performed at the low energy model scale Ql, and the 
resulting distributions are evoluted to the experimental large momentum scale in terms 
of the Altarelli-Parisi equation. Although next to leading order contributions of the QCD 
evolution are small pl| , it is very difficult to estimate uncertainties for the use of the QCD 



perturbation below IGeV. 

Our results are in reasonable agreements with experimental data, except for the large 
X region. In this region, the struck quark carries a large momentum > IGeV, and the 
NJL model is not designed for the momentum > ~ IGeV^. Generally, the phe- 
nomenological quark wave functions, e.g. Isgur-Karl model(Harmonic oscillator type) or 
the MIT bag model, do not include such high momentum components. We should improve 
the behavior of the structure function at the large x, by taking into account the quark 
correlation in the high momentum region. Higher twist contributions are also expected to 
change the shape of the distribution function at the large x. Indeed, experimental analysis 
of the scaling violation indicates that higher twists give non-negligible contributions to 
the nucleon structure functions in the large x region. If the pion is a collective qq state, 
twist-4 contributions may be larger than those of the nucleon. 
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The extraction of the flavor symmetry breaking effects in DIS is also important. By 
studying the differences between quark distributions of various flavors, we obtain the 
information about the flavor dependence of their internal wave functions in hadrons, e.g. 
the nucleon structure function gives the spin-flavor structure of the valence quarks for 
the u- and d-quarks||35|. The kaon structure function provides the valuable data on the 
strange sector. As we have discussed for the kaon, the valence strange quark may carry 
larger momentum fraction than the up or down quark in the kaon. Comparing the u-quark 
distribution in the kaon uk{x) with that in the pion Ut^{x), the NJL model calculation 
indicates the dominance of M7r(x) at the large x, and is consistent with the available data. 
This is due to the strong quark correlation in the pion, namely, the binding energy of the 
u-quark in the pion is larger than that in the kaon. 

The investigation of the semi-inclusive lepton nucleon scattering is very interesting as 
an alternative method to extract the virtual meson contributions to the nucleon structure. 
From them, we can obtain the off-shell dependence of the quark distribution function in the 
virtual pion. It also helps deeper understandings of the SU (2) flavor symmetry breaking 
in the nucleon sea. The calculated distribution function shows rather large momentum 
dependence. The pion structure function decreases as the spacelike momentum of the 
pion t becomes large, particularly for x > 0.4. We have computed the cross section with 
the use of the off-shell pion structure function and the on-shell one, and found the cross 
section of the off-shell t-dependent case is much smaller than the t = fixed case for 
X > 0.3. In order to compare our calculations with the forthcoming experiment, we ought 
to study this process more rigorously with the improvement of the large x behavior of 
the distribution function, e.g. estimation of contributions from other mesons and nucleon 
resonances, or inclusion of the interference of the vr and rj mesonfl^. Such a study is 



under investigation p6 



Here, we comment on the sea quark distributions in the NJL model. The lowest 
order diagram for the sea quark distributions is shown in Fig. 15 (a), corresponding to 
Fig.l (a), which is very important for the calculation of the quark self-energy and hence 
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the spontaneous chiral sjTiimetry breaking. As discussed by Bernard et a/. |]T^, however, 
this sea quark diagram does not contribute to the structure function in the NJL model, 
since it has no imaginary part. In order to get the 'meaningful' sea quark distribution 
in the model, we must calculate next to leading order diagram shown in Fig. 15 (b), 
which are usually neglected within the Hartree approximation. This inconsistency of the 
approximation is a defect of the NJL model. It can be shown that the contribution of the 
next to leading order diagram to the quark self-energy is about 10% of the total constituent 
mass. Therefore, the sea quark distribution is expected to share the momentum fraction 
of the same magnitude in the pion at the low energy model scale . 



Our approach in terms of the NJL model achieves qualitative agreements with the 
existing experimental data. The calculated quark distributions reflect the low energy 
quark dynamics of the model, i.e. the single-particle energy of the valence quark and 
the quark correlation in mesons. Such agreements may indicate the possible connection 
between the low energy quark model and the deep inelastic phenomena. 

We thank Prof. Y. Mizuno for valuable discussions about the semi-inclusive electron 
nucleon scattering. 
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Table 1 Meson Properties 
Parameters 

mu,d = 5.5MeV, rris = ISOMeV, A = 900MeV, 
GsA'^ = 2G'i/A2 = 19, a = O.lbGeV^ 





Theoretical values 


Empirical values 


Pion mass 


140MeV 


140MeV 


fn 


93MeV 


93MeV 


Kaon mass 


498MeV 


495MeV 


fK 


96MeV 


114MeV 


Mu 


370MeV 






548MeV 




<ijji/j > 


(-250Me1/)3 


- {-2bQMeVf 


9-Kqq 


3.92 
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Figure Captions 

Fig. 1 (a) : The Dyson equation for tlie quark propagator. Tlie tliick solid line represents 
the dressed constituent quark, and the thin solid line the current quark. 



Fig. 1 (b) : The Bethe-Salpeter equation for quark-antiquark scattering. T is the quark- 
antiquark T-matrix with the total momentum q. K, denotes the interaction kernel dis- 
cussed in the text. 

Fig. 2 : The forward scattering amplitude ("handbag diagram") of the pion. The solid 
line represents the quark, and dashed line the meson. The virtual photon is depicted by 
the wavy line. For details and notation, see text. 

Fig. 3 : The quark distributions of the pion at the low energy model scale = as a 
function of the Bjorken x. 

Fig. 4 : The quark distributions of the kaon at the low energy model scale = Ql- The 
u-quark distribution is depicted by the solid curve, and the s-quark distribution by the 
dashed curve. 

Fig. 5 : The valence quark distribution of the pion at = 20GeV^ (solid curve) as 
a function of the Bjorken x, in which we use the model scale Qq = 0.25GeV^. The 



experimental fit|2J] is depicted by the dashed curve. 



Fig. 6 : The valence quark distribution of the pion at = 20GeV^ with the model scale 
Ql = O.TGeV'^. The notations are the same as those in Fig. 5. 

Fig. 7 : The pion structure function at = 20GeV^. The experimental data are taken 
from the NA3 experiment ||12|| . The theoretical prediction of the NJL model is depicted 
by the solid curve. Here, we use F^{x) = Kxq{x) with the sea quark distributions being 
neglected, where we take the K-factor, K = 1.5 pi]. 
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Fig. 8 : The valence quark distributions in the kaon at = 20GeV'^ with the model 
scale Ql = 0.25GeV^. The u-quark distribution is depicted by the solid curve, and the 
s-quark distribution by the dashed curve. 



Fig. 9 : The ratio of kaon to pion valence u-quark distributions Uk{x) /ut^{x) at = 
20GeV^. The theoretical result is depicted by the solid curve. The closed circles with error 
bars are taken from the Drell-Yan experiment [T^]. The open circle line with the dotted 
area, which indicates the error, is obtained by taking a ratio of the parameterization in 



ref. where the pion and kaon structure functions are determined by the data on the 
large transverse momentum ttq productions. For comparison, the results with Ms/Mu =2 
and 1 are shown by the dashed and dotted curves, respectively. 

Fig. 10 : The valence quark distribution in the rho meson at the model scale as a function 
of the Bjorken x. 

Fig. 11 : The semi- inclusive lepton nucleon scattering. The solid line represents the 
electron and the doubled line the nucleon. The virtual pion and the virtual photon are 
depicted by the dashed and wavy hues. 

Fig. 12 : The off-shell dependence of the quark distribution in the pion. The results at 
t = ml (on-shell), and t = —0.1, —0.25, —O-bGeV"^ are shown. 

Fig. 13 : da/dp{e + p e' + n + X) as a function of p with El = SOGeV. In Fig. 13 (a), 
the triple Regge formula is used for the ttNN vertex, and the dipole cutoff in Fig. 13 (b). 
The theoretical result is depicted by the solid curve. The results using the distribution 
of Sutton et al. and Gliick et al. are also shown by the dashed and dotted curves. 



respectively. 

Fig. 14 : da/dx{e + p e' + n + X) as a function of the Bjorken x with El = 30GeV. 
The triple Regge formula is utilized in Fig. 14 (a), whereas the dipole formula in Fig. 14 
(b). The notations are the same as those in Fig. 13. The result in terms of the calculated 
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structure function with no off-shell dependence is also shown by the dash-dotted curve. 

Fig. 15(a) : The lowest order sea quark diagram in the NJL model. The solid line denotes 
the quark, and the wavy line the virtual photon. 

Fig. 15(b) : The next to leading order sea quark diagram in the NJL model. The dashed 
line indicates the pion. 
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